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total ionization cross sections for HJ, NJ, 
CD*, and CO* from the modified BEB 
formula are compared to available 
experimental data. The theory is in good 
agreement with experimental data on Hj, 
Nj, and CD*, but the experimental data 
by Belie et al. on NJ, and CO* are lower 
than the theory. The large difference 
between the theory and experiment on CO* 
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the dissociative ionization channel, 
CO*— >C*+0*, which was not included in 
the experiment. 
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1. Introduction 

Ionization cross sections for atomic and molecular 
ions are among the essential data needed in the study of 
plasmas in fusion devices as well as in a diverse class 
of mass spectrometric studies of molecules. The 
Binary-Encounter-Bethe (BEB) model [I] has been 
successful in generating reliable total ionization cross 
sections of small and large molecules [2-6]. The BEB 
model combines a modified form of the Mott cross 
section with the asymptotic form of the Bethe theory 
(i.e., high incident energy T) for electron-impact ioniza- 
tion of a neutral atom or molecule. 

In this article, we report a slight modification of the 
original BEB formula in order to adapt it to molecular 
ions. Almost all stable molecular ions are singly 
charged. The present modification has been found to 
reproduce known total ionization cross sections of 
singly charged ions, both atomic and molecular. We 
summarize the theory in Sec. 2, compare the theoretical 
cross sections to available experimental data in Sec. 3, 
and present our conclusion in Sec. 4. 



2. Summary of Theory 

The BEB formula for ionizing an electron from a 
molecular orbital of a neutral molecule by electron 
impact has a simple form [1]: 
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where t=TIB, u= U/B, S = 4'nai NR^IB^, a^ is the 
Bohr radius (= 0.5292 A), R is the rydberg energy 
(= 13.6057 eV), T is the incident electron energy, and 
N, B, and U are the electron occupation number, the 
binding energy, and the average kinetic energy of the 
orbital, respectively. 

In Eq. (1), the terms in the square brackets have 
sound theoretical origin based on the Mott theory and 
the Bethe theory. However, the denominator t+u+l is 
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based on a plausible, but less rigorous argument, i.e., the 
effective kinetic energy of the incident electron seen by 
the bound target electron should be the incident electron 
energy T plus the potential energy U+B of the target 
electron [7]. Hence the T in the denominator of the 
original Mott and Bethe theories was replaced by 
T+U+B, or t+u+l in Eq. (1), where B is used as the 
energy unit. This modification can be considered quali- 
tatively as a simple scheme to introduce T-dependent 
electron correlation between the incident and target 
electrons. 

The net effect of using t+u+l instead of t in the 
denominator of Eq. (1) is to reduce substantially the 
cross section near the ionization threshold. This modifi- 
cation was found not only to be effective but also 
absolutely necessary to bring the theory in agreement 
with reliable experimental ionization cross sections near 
the threshold for many neutral atoms and molecules. 

For singly charged ions, however, we noticed that this 
choice of the denominator reduced ionization cross 
sections too much near the threshold. This is qualita- 
tively understandable because the Bethe theory and the 
Rutherford cross section, which is the Mott cross section 
without the electron exchange terms, increase when the 
target is an ion. These cross sections increase as q^, 
where q is the net charge of the target [8]. There are, 
however, no quantitative guidelines to modify the BEB 
model for ionized targets. We found that when we 
replaced u+l by (u+l)/! in Eq. (1), the slightly 
modified BEB formula reproduced known ionization 
cross sections for singly charged atomic and molecular 
ions, i.e.. 
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Equation (2) for singly charged ions is as simple as the 
BEB formula for neutral targets, Eq. (1), and does not 
require any more input data than the original BEB 
formula. 

In addition, we also have a more detailed version of 
the BEB formula when a quantity Q related to the 
continuum oscillator strength, df/dE with photon 
energy E is known for each orbital [1]: 
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The value of Q is known for single orbital targets such 
as H, He, H-like ions, HJ, and H2 from either theory or 
photoionization experiments. We shall refer to this form 
of the cross section as the BEQ formula. The BEB 
formula, Eq. (1), was obtained by assuming Q=\ when 
df/dE for each orbital was not available. 

3. Results and Discussion 

We compare the BEB cross sections for HJ, N2, CD*, 
and CO* based on Eq. (2) to available experimental data 
[9-14] in Figs. 1-4. Most molecular constants B and U 
were calculated using the 6-3 1 1 -HG(d,p) basis set in the 
GAMESS code [15]. Since we were aware that electron 
correlation could increase the theoretical cross sections 
primarily through smaller B values [5], we obtained an 
alternative set of B and U for NJ orbitals with low 
binding energies from complete-active-space (CAS) 
wave functions using the MOLPRO code [16]. Also, 
our experience with the BEB model indicated that it 
was important to use an accurate value of the lowest 
ionization energy (IE) to obtain reliable cross sections 
near the threshold. For this purpose, we found the IE 
from the frozen-core coupled cluster theory to be suffi- 
cient [17-18] (single- and double-excitation operators 
[19] with a perturbative estimate of the contribution 
from the triple-excitation operator [20], using Dun- 
ning' s correlation-consistent valence-triple-zeta basis 
set [21]). Theoretical results and comparisons to exper- 
imental data presented in this article are also available 
on a NIST web page for electron-impact ionization 
cross sections at http://physics.nist.gov/ionxsec. 

In general, when an electron collides with a molecu- 
lar ion we get 



e +AB*^AH-B*H-e- or A^-nB-ne , 



e +AB*^AB* 
or AH-B**-i-2e^ 



■-i-2e-; A*H-B* + 2^-; A**-i-BH-2e- 



(«) 



(b) 



When Q is known, then Eq. (2) is replaced by 



Processes (a) are dissociation without ionization, 
while processes (b) are the ionizing events described by 
the BEB model. The model calculates the sum of all 
processes in (b) that lead to the ejection of a bound 
electron. The model can also account for double or 
higher multiple ionization resulting from the Auger de- 
cay of inner-shell holes [6], although the model is un- 
able to account for the simultaneous ejection of two or 
more electrons from a single molecular orbital. On the 
experimental side, it is difficult to distinguish processes 
(fl ) from processes ib) in a. single experiment unless a 
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Fig. 1. Comparison of the BEB and BEQ cross sections to experiment for HJ. Solid curve, 
BEQ cross section; dashed curve, BEB cross section; circles, experimental data by Peart and 
Bolder [9]. The ordinate is the ionization cross section, and the abscissa is the incident 
energy. The experimental uncertainties quoted by Peart and Bolder include 
uncertainties from both systematic and random effects. 
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Fig. 2. Comparison of the BEB cross section to experiments for Nj. Solid curve, BEB cross 
section using the constants from the restricted open-shell Hartree-Fock wave function; 
dashed curve, BEB cross section using the constants from the complete-active-space SCF 
wave functions; circles, experimental data by Peterson et al. [10]; and triangles, data by Belie 
et al. [12] and Siari et al. [11]. The two sets of experimental data represent slightly different 
physical quantities (see text). See Fig. 1 caption for other legend. 
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Fig. 3. Comparison of BEB cross section to experiment for CD*. Solid curve, BEB cross 
section for CH*; circles, experimental data by Djuric et al. [13] for CD*. See Fig. 1 caption 
for other legend. 
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Fig. 4. Comparison of BEB cross section to experiment for CO*. Solid curve, BEB cross 
section using a restricted open-shell Hartree-Fock wave function; dashed curve, BEB cross 
section using an unrestricted open-shell Hartree-Fock wave function; triangles, experimental 
data by Belie et al. [14]. See Fig. 1 caption for other legend. 



288 



Volume 105, Number 2, March-April 2000 

Journal of Research of the National Institute of Standards and Technology 



coincidence measurement of all products is performed. 
The usual experimental procedure is to measure the 
cross section for producing any ion, i.e., {a) + (b). 
Then, processes (a) are measured separately, and sub- 
tracted from the total ion production cross section. 
This subtraction introduces large uncertainties in the 
resulting experimental ionization cross sections. 

3.1 HI 

The molecular orbital constants B, U, and A^ for HI 
are listed in Table 1, and the BEB and BEQ cross 
sections are compared in Fig. 1 to the experimental data 
by Peart and Dolder [9]. The value of Q can be derived 
by subtracting the equivalent of Mfon for discrete excita- 
tions [22], which is known as M^xc and is equal to 
0.7081, from the equivalent of Ml,t, for total inelastic 
collisions [23], which is commonly known as ML 
and is equal to 0.8277. The resulting value of 
ML = Mi, -Ml, = 0. 1 196 leads to G = 0.5275 for H|. 
The comparison in Fig. 1 shows that the BEQ cross 
section is in slightly better agreement with experiment, 
as the BEQ formula contains more accurate information 
about the ionization continuum than the BEB formula. 

3.2 N| 

The molecular orbital constants B, U, and N for N| are 
listed in Table 1, and the BEB model is compared in 
Fig. 2 to the experimental data by Peterson et. al [10] 
and by Siari et al. [11] and Belie et al. [12]. The exper- 
imental data by Peterson et al. are the difference 
between the total cross section for the production of N* 
measured by Van Zyl and Dunn [24] and the dissociative 
excitation 



NlH-e-^N^H-N-ne- 



(6) 



measured by Peterson et al. Although this difference, 
which is shown in Fig. 2, is in excellent agreement with 
the BEB cross section, the experimental difference 
counts the symmetric dissociation 



N|-He-^N*-i-N*H-2e- 



(7) 



twice because this process produces two N* ions. How- 
ever, the experimental data by Belie et al. [12] in Fig. 2 
counts the symmetric ionization process [Eq. (7)] only 
once. The peak value of the cross section for the process 
[Eq. (7)] measured by Belie et al. [12] is about 0.25 Al 
If we reduce the data by Peterson et al. [10] by this 
amount, the lower half of their error bars will overlap 
with the data by Belie et al., which are estimated to have 
(5 to 6) % combined relative uncertainties due to sys- 
tematic and random effects [11,12]. 

The BEB model counts each ionizing event only once. 
However, the theory assumes that every energy transfer 
exceeding the lowest IE will result in ionization, as all 



Table 1. Molecular orbitals, electron binding energy B, kinetic 
energy U, electron occupation number N, and the dipole quantity 
Q for H2, NJ, CD*, and CO*. Wave functions: RHF stands for 
restricted Hartree-Fock, CAS means complete active space, and UHF 
is unrestricted Hartee-Fock. 
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" Ionization energies deduced from experimental data. 
Ionization energies from the coupled-cluster method [17-21]. 

binary-encounter type theories do. For molecules, this 
assumption is a poor one because molecules can dis- 
sociate without ejecting an electron as in the dissociative 
excitation process [Eq. (6)] even when it receives an 
energy transfer exceeding the lowest IE. In short, the 
BEB model is likely to overestimate the total ionization 
cross section. Our experience with CF4 [5] indicates that 
the BEB cross section derived from correlated molecular 
wave functions is likely to be (5 to 10) % higher at the 
peak than the correct peak value. 

3.3 CD* 

The molecular orbital constants for CD* are listed in 
Table 1, and the BEB cross section is compared in 
Fig. 3 to the experimental data by Djuric et al. [13]. The 
experimental ionization cross section was deduced by 
subtracting the cross section for dissociative excitation 
measured by Forck [25] from the total D* production 
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cross section measured by Djuric et al. [13]. The BEB 
cross section is insensitive to isotopes, and the theoret- 
ical cross section is for CH*, while the experiment was 
done with deuterium. No uncertainty estimates are 
available for the experimental data. 

3.4 CO* 

The molecular orbital constants for CO* are listed in 
Table 1, and the BEB cross section is compared in 
Fig. 4 to the experimental data by Belie et al. [14]. The 
molecular constants for the BEB cross section were 
determined by using two different molecular wave 
functions for the molecular ion. Both sets of constants 
produce basically the same cross section. The experi- 
ment by Belie et al. [14] left out an important ionization 
channel: 

C0*H-e-^C*-i-0*-i-2e- . (8) 

Their experiment measured the cross section for the 
production of doubly charged ions, CO**, C**, and O**. 
As is shown in Fig. 4, the experimental cross section by 
Belie et al. [14], who estimate the combined relative 
uncertainties from systematic and random effects to be 
(5 to 8) %, is almost an order of magnitude lower than 
the present theory. The neutral molecules CO and N2 
are isoeleetronic, and their ionization cross sections 
are similar in magnitude and shape [2]. Hence, we 
expect ionization cross sections of CO* and N2 to be 
also similar in shape and magnitude. Such similarity is 
observed in the BEB cross sections in Figs. 2 and 4. The 
discrepancy between the theory and experiment in 
Fig. 4 indicates that the missing ionization process 
[Eq. (8)], which is equivalent to the symmetric dissocia- 
tive ionization process [Eq. (7)] in a homopolar di- 
atomic molecular ion, is the dominant mode for ioniz- 
ing CO*. 

4. Conclusion 

We have shown that a simple modification, Eq. (2), 
of the BEB formula for neutral molecules, Eq. (1), 
produces reliable electron-impact total ionization cross 
sections for HJ and CD*. On the other hand, the exper- 
imental data on NJ by Belie et al. [12] are lower than 
the modified BEB cross section, while the experimental 
data deduced from the experiments by Van Zyl and 
Dunn [24] and by Peterson et al. [10] are in good agree- 
ment with the present theory taking into consideration 
the fact that the data by Peterson et al. counted the 
symmetric ionization process [Eq. (7)] twice, whereas 
the BEB cross section with the CAS wave function is 
likely to be (5 to 10) % too high as explained in Sec. 3.2. 

Belie et al. [14] measured only part of the ionization 
cross section for CO* that produces doubly charged 
ions. Judging from the large difference (see Fig. 4) 
between the present theory and the data by Belie et al. 
[14], it is likely that the major part of the ionization 



cross section of CO* is in the "symmetric" dissociative 
ionization process [Eq. (8)], unlike the case for N2, 
whose cross section for nonsymmetrie ionization is 
about twice higher than that for the symmetric ioniza- 
tion [12]. 
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